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Abstract. We set new limits on the mass of the Milky Way, making use of the latest kinematic information for 
Galactic satellites and halo objects. Our sample consists of 11 satellite galaxies, 137 globular clusters, and 413 
field horizontal-branch (FHB) stars up to distances of 10 kpc from the Sun. Roughly half of the objects in this 
sample have measured proper motions, permitting the use of their full space motions in our analysis. In order 
to bind these sample objects to the Galaxy, their rest-frame velocities must be lower than their escape velocities 
at their estimated distances. This constraint enables us to show that the mass estimate of the Galaxy is largely 
affected by several high- velocity objects (Leo I, Pal 3, Draco, and a few FHB stars), not by a single object alone 
(such as Leo I), as has often been the case in past analyses. We also find that a gravitational potential that gives 
rise to a declining rotation curve is insufficient to bind many of our sample objects to the Galaxy; a possible 
lower limit on the mass of the Galaxy is about 2.2 x 10 12 M @ . To be more quantitative, we adopt a Bayesian 
likelihood approach to reproduce the observed distribution of the current positions and motions of the sample, 
in a prescribed Galactic potential that yields a flat rotation curve. This method enables a search for the most 
likely total mass of the Galaxy, without undue influence in the final result arising from the presence or absence 
of Leo I, provided that both radial velocities and proper motions are used. Although the best mass estimate 
depends somewhat on the model assumptions, such as the unknown prior probabilities for the model parameters, 
the resultant systematic change in the mass estimate is confined to a relatively narrow range of a few times 
lO n M , owing to our consideration of many FHB stars. The most likely total mass derived from this method 
is 2.5+";o x 10 12 M© (including Leo I), and 1.8+Q 7 x 10 12 M Q (excluding Leo I). The derived mass estimate of 
the Galaxy within the distance to the Large Magellanic Cloud (~ 50 kpc) is essentially independent of the model 
parameters, yielding 5.5lg 2 x Mq (including Leo I) and 5.4^0 4 x 10 11 Mq (excluding Leo I). Implications 
for the origin of halo microlensing events (e.g., the possibility of brown dwarfs as the origin of the microlensing 
events toward the LMC, may be excluded by our lower mass limit) and prospects for more accurate estimates of 
the total mass are also discussed. 
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ies, and the observed gravitational lensing of background 
sources (e.g., Binney & Tremaine 1987). A determina- 
tion of the extent over which such dark-matter-dominated 
mass distributions apply for most galaxies, including our 
own, is of great importance for understanding the role of 
dark matter in galaxy formation and dynamical evolution. 
In particular, the mass estimate of the Galaxy is closely 
relevant to understanding the origin of the microlensing 
events toward the Large Magellanic Cloud (LMC) (e.g., 
Alcock et al. 2000; Alcock et al. 2001). 

While mass estimates of external galaxies can, in prin- 
ciple, be obtained in a straightforward fashion using var- 
ious dynamical probes, the total mass of the Galaxy re- 



1. INTRODUCTION 

Over the past decades, various lines of evidence have re- 
vealed that the mass density in the Milky Way is largely 
dominated by unseen dark matter, from the solar neigh- 
borhood to the outer reaches of the halo (e.g., Fich & 
Tremaine 1991). Moreover, the presence of a dark compo- 
nent similar to that found in our own Galaxy appears to be 
a generic feature in external galaxies, as inferred from, e.g., 
flat rotation curves in their outer parts, the presence of 
(a gravitationally bound) hot plasma in early-type galax- 
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mains rather uncertain, primarily due to the lack of ac- 
curate observational information for tracers located in its 
outer regions, where the dark matter dominates. The pre- 
cise shape of the outer rotation curve, as deduced from 
H II regions and/or H I gas clouds (e.g., Honma & Sofue 
1997), is still uncertain because its determination requires 
knowledge of accurate distances to these tracers (Fich & 
Tremaine 1991). Also, interstellar gas can be traced only 
up to ~ 20 kpc from the Galactic Center, and hence pro- 
vides no information concerning the large amount of dark 
matter beyond this distance. 

The most suitable tracers for determination of the 
mass distribution in the outer halo of the Galaxy are 
the distant luminous objects, such as satellite galaxies, 
globular clusters, and halo stars on orbits that explore 
its farthest reaches (e.g., Miyamoto et al. 1980; Little 
& Tremaine 1987; Zaritsky et al. 1989; Kochanek 1996; 
Wilkinson & Evans 1999, hereafter WE99). However, the 
limited amount of data presently available on the full 
space motions of these tracers, and the small size of the 
available samples, have stymied their use for an accurate 
determination of the Galaxy's mass. In particular, most 
previous mass estimates (except for WE99, see below) de- 
pend quite sensitively on whether or not a distant satellite, 
Leo I, is bound to the Galaxy. Leo I has one of the largest 
radial velocities of the known satellites, despite its being 
the second most distant satellite from the Galaxy (Mateo 
1998; Held et al. 2001). As a consequence, estimates of the 
total mass of the Galaxy are much more uncertain (by as 
much as an order of magnitude) than, for instance, the 
value of the circular speed in the solar neighborhood (Kerr 
& Lynden-Bell 1986; Fich & Tremaine 1991; Miyamoto & 
Zhu 1998; Mcndez et al. 1999). 

Recently, by making use of both the observed radial ve- 
locities and proper motions of six distant objects, WE99 
demonstrated that the use of full space motions can pro- 
vide a reliable mass estimate of the Galaxy without be- 
ing largely affected by the presence or absence of Leo I. 
They also argued that the primary uncertainties in their 
mass estimate arose from the small size of the data set 
and the measurement errors in the full space motions, es- 
pecially the proper motions. This work motivated us to 
investigate a much larger data set, with more accurate 
kinematic information, to set tighter limits on the mass 
of the Galaxy. Specifically, as we show below, there are 
two objects among the WE99 sample (Draco and Pal 3) 
that have relatively large velocity errors, yet still play cru- 
cial roles in a determination of the Galaxy's mass, so the 
addition of more (and better) data is important. 

Over the past few years, the number of distant satellite 
galaxies and globular clusters with available proper mo- 
tions has gradually increased (e.g., Mateo 1998; Dinescu 
et al. 1999; Dinescu et al. 2000; Dinescu et al. 2001). In 
addition, another tracer population that is suitable for 
exploring mass estimates of the Galaxy has become avail- 
able from the extensive compilation of A-type metal-poor 
stars by Wilhelm et al. (1999b), which provided radial 
velocity measurements, as well as estimates of the phys- 



ical parameters of these stars (e.g., [Fe/H], T e fj, log g). 
Among the Wilhelm et al. sample, the luminous FHB 
stars are the most useful mass tracers, both because of 
their intrinsic brightness, and the fact that accurate dis- 
tance determinations can be inferred from their absolute 
magnitudes on the horizontal branch (e.g., Carretta et al. 
2000). Moreover, there exist proper- motion measurements 
for many of these stars, provided by both ground- and 
space-based proper-motion catalogs (Klemola et al. 1994; 
Roser 1996; Platais et al. 1998; Hog et al. 2000), from 
which full space motions may be derived. 

In this paper we re-visit the mass determination of 
the Galaxy, based on a sample of 11 satellite galaxies, 
137 globular clusters, and 413 FHB stars, out of which 5 
satellite galaxies, 41 globular clusters, and 211 FHB stars 
have measured proper motions. First, we investigate the 
lowest possible mass that the Galaxy may have, adopting 
the requirement that the rest-frame velocities of observed 
sample objects be less than their escape velocities at their 
present distance from the Galactic center (e.g., Miyamoto 
et al. 1980; Carney et al. 1988). Secondly, the most likely 
mass of the Galaxy is calculated, based on a Bayesian like- 
lihood analysis that seeks to reproduce both the current 
positions and velocities of the sample objects (e.g., Little 
& Tremaine 1987; Kochanek 1996; WE99). Because our 
sample of tracers is, by far, the largest and most accu- 
rate one presently available, it is possible to place more 
reliable limits on the total mass of the Galaxy. In Sect. 2 
we describe our sample objects and the assembly of their 
kinematic data. In Sect. 3 we discuss the influence of the 
adopted membership of our sample on the mass determi- 
nation. In Sect. 4 we obtain the most likely total mass 
of the Galaxy based on a Bayesian likelihood analysis. In 
Sect. 5 we discuss implications for the origin of the halo 
microlensing events toward the LMC and the mass esti- 
mate of the Local Group, and consider the prospects for 
more obtaining more accurate estimates of the total mass 
of the Galaxy in the near future. 



2. DATA 

We consider a sample of objects that serve as tracers of 
the Galactic mass distribution consisting of 11 satellite 
galaxies, 137 globular clusters, and 413 FHB stars. In the 
case of the satellite galaxies, all of the basic information 
for our kinematic analysis, i.e., positions, heliocentric dis- 
tances, and heliocentric radial velocities, are taken from 
the compilation of Mateo (1998). For the globular clus- 
ters, we adopt the information provided by Harris (1996), 
including their positions and heliocentric radial velocities, 
their metal abundances, [Fe/H], and the apparent magni- 
tude of the clusters' horizontal branch (HB). The catalog 
of Wilhelm et al. (1999b) is our source of similar informa- 
tion for the FHB stars. We obtain an internally consistent 
set of distance estimates for the globular clusters and the 
FHB stars from the recently derived relationship between 
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the absolute magnitude of the HB, My (HB), and [Fe/H], 
by Carretta et al. (2000), 

My (HB) = (0.18 ±0.09) ([Fe/H] + 1.5) + (0.63 ±0.07) .(1) 

Clearly, we have assumed that there is no large offset be- 
tween the absolute magnitudes of FHB stars and their 
counterpart HB stars in the globular clusters (a view also 
supported by the recent work of Carretta et al. 2000). Fig. 
1 shows the spatial distribution of the globular clusters, 
satellite galaxies, and FHB stars on the plane perpen- 
dicular to the Galactic disk, where X axis connects the 
Galactic center (X=Q) and the Sun (X=8.0 kpc). The 
filled and open symbols denote the objects with and with- 
out proper-motion measurements, respectively. Satellite 
galaxies are the most distant tracers, with Galactocentric 
distances, r, greater than 50 kpc. The globular clusters 
extend out to almost r — 40 kpc, while the present sam- 
ple of FHB stars are confined to locations within 10 kpc 
of the Sun. Thus, our sample objects are widely, though 
not uniformly, distributed throughout the volume of the 
Galaxy. 

Among these sample objects there exist proper-motion 
measurements for 5 of the satellite galaxies, 41 of the glob- 
ular clusters, and 211 of the FHB stars. The proper motion 
data for the LMC, Sculptor, and Ursa Minor are taken 
from WE99, whereas those for Sagittarius and Draco are 
taken from Irwin et al. (1996) and Scholz & Irwin (1994), 
respectively. The proper motions for most of the globu- 
lar clusters have been compiled by Dinescu et al. (1999). 
We adopt the data from this source, except for two globu- 
lar clusters with recently revised proper-motion measure- 
ments (NGC 6254: Chen et al. 2000; NGC 4147: Wang et 
al. 2000), and for three additional globular clusters com- 
piled recently (Pal 13: Siegel et al. 2000; Pal 12: Dinescu et 
al. 2000; NGC 7006: Dinescu et al. 2001). Proper motions 
for 211 of the FHB stars in the Wilhelm et al. (1999b) 
sample are available from one or more existing proper- 
motion catalogs. These include the STARNET Catalog 
(Roser 1996), the Yale-San Juan Southern Proper Motion 
Catalog (SPM 2.0: Platais et al. 1998), the Lick Northern 
Proper Motion Catalog (NPM1: Klemola et al. 1994), and 
the TYCHO-2 Catalog (Hog et al. 2000). Many of these 
FHB stars have been independently measured in two or 
more catalogs, so that by combining all measurements one 
can reduce the statistical errors, as well as minimize any 
small remaining systematic errors in the individual cat- 
alogs, as was argued in Martin & Morrison (1998) and 
Beers et al. (2000). 

We estimate average proper motions, < fi >, and their 
errors, < cr M >, weighted by the inverse variances, as 
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Fig. 1. Spatial distributions of satellite galaxies (squares), 
globular clusters (circles), and FHB stars (triangles) on 
the plane perpendicular to the Galactic disk, where the 
X axis connects the Galactic center {X =0) and the Sun 
(X=8.0 kpc). The filled and open symbols denote the ob- 
jects with and without available proper motions, respec- 
tively. The plus sign in panel (b) denotes the position of 
the Sun, (X,Y) = (8.0,0). 



< M > = £>/<)/(£ V<), 

i=l i=l 

<*„> = or i/O- 1 / 2 , 



»=i 



( ) where n denotes the number of measurements. Table 1 
(available at the CDS) lists these compilations, as well 

^ as the estimated distances to the FHB stars, where r 
and RV denote the Galactocentric distances and hclioccn- 
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trie radial velocities, respectively. Typical errors in the 
reported proper-motion measurements range from 1 to 
~ 5 mas yr" 1 for individual field stars, whereas those for 
the satellite galaxies and globular clusters are about 0.3 
mas yr _1 and 1 mas yr _1 , respectively. 

We assume a circular rotation speed for the Galaxy 
of Vlsr = 220 km s -1 at the location of the Sun (i.e. 
Rq = 8.0 kpc along the disk plane) and a solar motion of 
(U, V, W)= (-9, 12, 7) km s" 1 (Mihalas & Binney 1981), 
where U is directed outward from the Galactic Center, V 
is positive in the direction of Galactic rotation, and W is 
positive toward the North Galactic PoleQ. We then calcu- 
late the space motions of our tracers, as well as the errors 
on these motions, fully taking into account the reported 
measurement errors in the radial velocities of the individ- 
ual satellite galaxies (typically a few km s _1 ), adopting a 
typical radial-velocity error for other objects (10 km s" 1 ), 
the measurement errors assigned to the proper motions 
of each object (when available, adopting a mean error for 
the source catalog when not), and distance errors for the 
satellite galaxies (10 % relative to the measured ones), or 
as obtained from Eq. 1 for the globular clusters and FHB 
stars. 

It is worth noting that the reported proper motions of 
the FHB stars in our sample may yet contain unknown 
systematics with respect to their absolute motions in a 
proper reference frame; this caution applies to the globu- 
lar clusters and satellite galaxies as well. It is an important 
goal to make efforts to reduce the systematic, as well as 
random, errors in the proper motions upon which studies 
of Galactic structure and kinematic studies are based, us- 
ing much higher precision astrometric observations than 
have been obtained to date. 
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Fig. 2. Rotation curves for Model A and Model B, param- 
etcrizations of the mass distributions considered in this 
paper. See the text for more information on the nature of 
these models. 



of elucidating the role of sample selection in setting limits 
on the mass of the Galaxy. 

To this end we consider two different mass models, 
in order to investigate the difference in the mass limits 
obtained by the use of different potentials. Our models, 
hereafter referred to as Model A and B, are the same as 
those adopted in WE99 and Johnston et al. (1995) (and 
also used by Dinescu et al. 1999), respectively. 

Model A has spherical symmetry, and results in a flat 
rotation curve in the inner regions of the Galaxy. The 
gravitational potential and mass density are given as: 



3. ROLES OF THE SAMPLE IN THE MASS 
LIMITS 

If we model the Galaxy as an isolated, stationary mass 
distribution, and assume that all of our tracer objects are 
gravitationally bound to it, then the rest-frame velocities 
of all objects, Vrf, must be less than their escape ve- 
locities, V csc = y/Ztfi, where ip denotes the gravitational 
potential of the Galaxy. A number of previous researchers 
have adopted this method for estimation of the mass of the 
Galaxy (e.g., Fricke 1949; Miyamoto et al. 1980; Carney 
et al. 1988; Leonard & Tremaine 1990; Daupholc & Colin 
1995). Note that the mass estimate obtained in this fash- 
ion is largely dominated by only a small number of high- 
velocity objects, hence the mass that is derived depends 
rather sensitively on the selection criteria adopted for such 
objects. In the present section, instead of deriving an exact 
mass estimate, we follow this procedure for the purpose 



1 Dehnen & Binney (1998) derived the solar motion of 
(U, V, W) = (-10.0, 5.3, 7.2) km s" 1 based on Hipparcos data. 
Because the difference between this and the currently adopted 
solar motion is only a few km s _1 , it gives little influence on 
the mass estimate. 
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(4) 



where a is the scale length of the mass distribution, and 
M is the total mass of the system. The central density 
of this model is cusped (like r~ 2 ) and falls off as r~ 5 for 
r 3> a. Using Eq. (4), the circular rotation speed is given 
as V? = GM/ (r 2 + a 2 ) 1/2 , so by setting V c at r = R Q as 
Vlsr = 220 km s" 1 in our standard case, it follows that 
this model contains one free parameter, o, to obtain M. 

Model B consists of more realistic axisymmetric po- 
tentials with three components (the bulge, disk, and dark 
halo) that reproduce the shape of the Galactic rotation 
curve (Johnston et al. 1995). The bulge and disk compo- 
nents we adopt are those represented by Hernquist (1990) 
and Miyamoto & Nagai (1975) potentials, respectively. All 
of the parameters included in these potentials are taken 
from Dinescu et al. (1999) (see their Table 3). In order to 
obtain a finite total mass we assume the following modified 
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logarithmic potential (corresponding to an isothermal-like 
density distribution) for the dark halo component: 



r^log[l + (r/d) 2 ]-Vo, 
^ hal ° l ' ~ 1 -2wg^T^, at r>r cut 



at r < r cut 



p(r) = 



2v% 



3 + r/d 



4nGd 2 (1 + r/df ' 
where ipo is defined as 
^o = «gpog(l + c) + 2c/(l + c)] ) 



c = (r cut /dy 



(5) 



(6) 



(7) 



and we adopt v = 128 km s _1 and d = 12 kpc (Dinescu et 
al. 1999). This model contains one free parameter, namely 
the cutoff radius of the dark halo, r cut . Fig. 2 shows the 
rotation curves for < R < 20 kpc, provided by Model 
A with a — 200 kpc (thick solid line) and Model B with 
?"cut = 170 kpc (thin solid line), where both curves shown 
at R < 20 kpc remain unchanged as long as a, r cut 3> 20 
kpc. The circular speed at R — Rq is 220 km s _1 for both 
mass models. Also shown is the declining rotation curve 
with increasing radius, as obtained from Model A with 
a = 20 kpc and Vlsr = 211 km s _1 (dashed line). 

We plot, in Fig. 3a and 3b, the relationship between 
the derived escape velocities, V eS c, and the rest-frame ve- 
locities, Vrf, when we adopt Model A with a = 195 kpc 
and Model B with r cut = 295 kpc, respectively. For the 
objects without available proper motions (open symbols), 
we adopt the radial velocities alone as measures of Vrf, 
hence their estimated space velocities are only lower limits. 
The solid line denotes the boundary between the objects 
that are bound (below the line) and unbound (above the 
line) to the Galaxy, respectively, where the positions of the 
data points relative to the solid line depend on the choice 
of a or r cu t . As these figures clearly indicate, the selection 
of the smallest value of a or r cu t that places the sample 
inside the bound region, or equivalently the lower limit 
on the mass of the Galaxy, is controlled by a few high- 
velocity objects located near the boundary line at each 
respective radius (or corresponding ip). These objects in- 
clude Leo I (for which only radial velocity information is 
available), Draco, Pal 3, and several high- velocity FHB 
stars. This highlights the following important properties 
of the derived mass limits: (1) If the proper motions of 
all objects are unavailable, then the mass estimate sen- 
sitively depends on the presence or absence of Leo I, as 
has been noted in previous studies. (2) Compared to case 
(1), if the available proper motions of the satellite galax- 
ies and globular clusters are taken into account, the con- 
straint provided by Draco and Pal 3 is basically the same 
as that provided by Leo I. This may explain the result 
of WE99, which showed that the mass determination is 
made insensitive to Leo I if the proper motion data of 
satellite galaxies and globular clusters are taken into ac- 
count. However, as Fig. 3 indicates, the velocity errors for 
Draco and Pal 3 are quite large, so these objects place 
only weak constraints on the mass estimate. (3) If we con- 
sider the proper motions of FHB stars, then some FHB 
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Fig. 3. (a) The relation between escape velocities, V csc , 
and space velocities, Vrf, for Model A with a = 195 kpc 
and Vlsr = 220 km s _1 . The symbols are the same as 
those in Fig. 1. The solid line denotes the boundary be- 
tween the gravitationally bound and unbound objects - 
those in the region below the line are bound to the Galaxy. 
For the sake of clarity, velocity errors are plotted for only 
the high-velocity objects relevant to the mass estimate, 
(b) Same as panel (a) but for Model B with r cut = 295 
kpc. 



stars having high velocities provide the basically the same 
constraint on the Galaxy's mass as do Leo I, Draco, and 
Pal 3. These properties suggest that the inclusion of FHB 
stars with available proper motions is crucial, and that 
they provide constraints on the mass limit of the Galaxy 
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that depend on neither the inclusion or absence of Leo I 
nor on the large velocity errors for Draco and Pal 3. 

As mentioned above, a determination of the lower mass 
limit for the Galaxy, using escape velocities, inevitably de- 
pends on the selection of a few apparently high- velocity ob- 
jects from a much larger sample of tracers. While a mass 
estimate independent of this selection effect will be ob- 
tained in Sect. 4, we seek first to obtain a rough measure 
for the lower mass limit, i.e., the smallest a or r cut that 
encloses Pal 3, Draco, and the four highest-velocity FHB 
stars (for which we mark asterisks after their names in 
Table 1) inside the bound region, based on a weighted 
least-squares fitting procedure (weights being inversely 
proportional to the velocity errors). This exercise yields 
a = 195±gf kpc for Model A, and r cut = 295±f|| kpc for 
Model B. Using these values, the lower limits to the total 
mass, M, of the Galaxy may be given as 2.2t 1 'o x 1O 12 M 
for Model A and 2.2±f;f x 1O 12 M for Model B, respec- 
tively. Thus, the difference in the derived mass limits is 
not significant, as long as the rotation curve at outer 
radii is approximately constant at the adopted value of 
220 km s _1 . It also suggests that the flattened nature of 
the Model B potential, due to the presence of the disk 
component, does not affect the results significantly - the 
high-velocity tracers are located at large Galactocentric 
distances and/or their orbits largely deviate from the disk 
plane. 

In addition to the above experiments, we also consider 
a mass model that yields a declining rotation curve at 
outer radii, as was proposed by Honma & Sofuc (1997) 
from their H I observations. We adopt Model A with a = 
20 kpc and Vlsr = 211 km s _1 at R = R@, so as to 
yield the declining rotation curve at R > Rq (dashed line 
in Fig. 2), which is reminiscent of the result in Honma 
& Sofue (1997). Fig. 4 shows the Vrf vs. V csc relationship 
that follows from adoption of this model. As is evident, the 
total mass obtained from a model that leads to a declining 
rotation curve is insufficient to bind many of our sample 
objects to the Galaxy. 



4. MASS DETERMINATION BASED ON A 
BAYESIAN LIKELIHOOD METHOD 

4.1. Method 

To obtain a more quantitative measure of the mass of the 
Galaxy, we examine an alternative method that takes into 
account all of the positional and kinematic information of 
the sample objects, in contrast to the use of the high- 
velocity tracers alone, as in the previous section. In this 
approach, a phase-space distribution function of tracers, 
F, is prescribed for a specifically chosen ip, and the model 
parameters included in F and ip are derived so as to re- 
produce the presently observed positions and velocities of 
the tracers in the (statistically) most significant manner. 
The optimal deduced parameters relevant to then allow 
us to estimate the total mass of the Galaxy. This method 
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Fig. 4. The relation between escape velocities, V^ sc , and 
space velocities, Vrf, for Model A with a = 20 kpc and 
Vlsr = 211 km s _1 . In this case, the rotation curve de- 
clines with increasing radii, as shown in Fig. 2 (dashed 
line). Note that, if this situation were to apply, many of 
the sample objects would be unbound to the Galaxy. 

was originally proposed by Little & Tremaine (1987), and 
further developed by Kochanek (1996) and WE99. 

Based on the results presented in the previous section, 
we take Model A with spherical symmetry as the mass 
distribution of the Galaxy, which is sufficient for the fol- 
lowing analysis. For the sake of simplicity, and also for 
ease of comparison with the previous studies by Kochanek 
(1996) and WE99, the phase-space distribution function is 
taken to have the same anisotropic form as that adopted 
in these studies. That is, it depends on the binding energy 
per unit mass, e (= ip — v 2 /2), and the angular momentum 
per unit mass, I, in the following way: 



F{e,l)=l- 20 f{e) 



(8) 



where 



2/3-3/2 



7r 3/2 r [ m _ 1/2 + p]T[l - (3} de 

u ' V s I- „/.\/3-3/2+m 



I 



dip- 



dip r ' 



'-{e-0 



(9) 



where p s is the tracer density distribution, L is the gamma 
function, and m is an integer whose value is chosen such 
that the integral in Eq. (9) converges (e.g., Dejonghe 1986; 
Kochanek 1996). In the spherical model, this form of the 
distribution function yields equal velocity dispersions in 
the orthogonal angular directions, < >—< >, and a 
constant anisotropy (3=1— < v@ > / < v% > everywhere 
in the Galaxy. Our choice of to = 2 in Eq. (9) (to be in 
accord with the WE99 work) limits the allowed range for 
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Fig. 5. Cumulative number distribution, N(< r), of the 
distances of globular clusters and satellite galaxies (solid 
histogram) in comparison with model distributions (con- 
tinuous dashed and solid lines). See the text for additional 
information. 

the velocity anisotropy to— 1.5</3<1 when proper mo- 
tion data are considered, while the use of radial velocities 
alone sets no limit for tangential anisotropy [-co, 1]. 

For p s , we consider WE99's two models: (a) Shadow 
tracers following the mass density distribution obtained 
from Model A (Eq. 4), and (b) a power-law distribution 
as a function of r. The shadow-tracer model is given as: 



Ps(r) oc 



r 2 (r 2 + a 2 ) 3 / 2 



(10) 



where a s is the scale length. The power-law model with 
index 7 is given as: 



Ps(r) oc 4: • 
r ~ 



(11) 



It should be noted that, since shadow tracers may be trun- 
cated at the distance below the scale length of the mass 
distribution, the scale length of the tracers, a s , is gener- 
ally different from the scale length of the Galaxy's mass, 
a. 

Using the 27 objects (satellite galaxies and globular 
clusters) at r > 20 kpc, WE99 derived a s = 100 kpc and 
7 = 3.4 as the best fitting parameters for their spatial 
distribution. We re-examine a s and 7 using our sample of 
all satellite galaxies and globular clusters. Note that the 
FHB stars are excluded in this determination of a s and 
7, as they have not (yet) been completely surveyed over 
the Galactic volume. We obtain a s = 10 kpc and 7 = 3.3 
as the best-fit values, based on a simple K-S test of the 
observed vs. predicted distribution functions (see Fig. 5). 
If we exclude the globular clusters at r < 10 kpc, for 
which the spherical symmetry assumption may be ques- 
tionable due to the presence of the disk globular clusters, 
we obtain a s = 50 kpc and 7 = 3.4. Thus, a s depends 
sensitively on the adopted range of radius (i.e., on the 



selection of the sample), whereas 7 basically remains un- 
changed. Therefore, we focus our attention on the results 
obtained from the power-law representation for the tracer 
population; the shadow-tracer population is also examined 
for the purpose of comparison with WE99. To see the de- 
pendence of the mass estimate on these parameters, we 
obtain estimates for two values of 7 (3.4 and 4.0) and a s 
(100 kpc and the scale length of the mass distribution, a), 
respectively. We note that the FHB stars are also expected 
to follow a power-law form with 7 ~ 3.4, as inferred from 
other halo field stars (e.g., Preston et al. 1991; Chiba & 
Beers 2001). 

We calculate the likelihood of a particular set of model 
parameters (the scale length of the mass distribution, a, 
and the anisotropy parameter, 0) given the positions, rj, 



and radial velocities, 



or space velocities, 



using 



Bayes' theorem. The probability that the model param- 
eters take the values a and 0, given the data (r"j, i>( r )j) 
and prior information /, is: 

1 N 

P{a,0\n,v {r)i ,I) = -P{a)P{0)J{P{r i ,v (r)i \a,0) , (12) 



where N is the normalization factor (Kochanck 1996; 
WE99). The probabilities P(a) and P{0) denote the prior 
probability distributions in a and 0, respectively. Here, 
P(ri, V( r )i\a, 0) corresponds to the probability of finding 
an object at position rj moving with radial velocity W( r ),, 
or space velocity Vi , for a particular set of model parame- 
ters a and 0. The complete expressions for P(r"j, vi r \j\a, 0) 
are shown in Table 1 of WE99. To calculate this probabil- 
ity for the objects with full space velocities, we take into 
account their large errors relative to radial velocities alone 
(due to the observed proper-motion and assumed distance 
errors), by multiplying by an error convolution function of 
the form: 

P(n,Vi\a,(3) = J J dv a dv s Ei(v a )Ei(vg) 

x P(ri,Vi i0bs (v a ,v s )\a, 0) , (13) 

where (v a , vg) are the tangential velocities along the right 
ascension and declination coordinates, respectively, and 
Ei is the Lorentzian error convolution function, defined 
as: 



Ei{v) 



1 



2a\ 



V27R71 2(T 2 + (V- V obs f 



(14) 



where 01 is defined as <Ji — 0A77a for the calibrated error 
estimate a (see WE99). 

The prior probability in the velocity anisotropy, 0, is 
taken to be of the form P(0) oc 1/(3 - 2/3)", where n = 
and 2 correspond to a uniform prior and a uniform energy 
prior, respectively (Kochanek 1996; WE99). Larger values 
of n yield a larger weight towards radial anisotropy. For 
the prior probability in a, P(a), we adopt 1/a and 1/a 2 
(WE99). 
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Using the routine AMOEBA in Numerical Recipes 
(Press et al. 1992), we search for a set of model parameters, 
a and (3, that maximize the probability P(a, /3|rj, U( r )j, /). 
The total mass of the Galaxy, M, is then derived from the 
parameter a. 



4.2. Results 

Initially, we apply the Bayesian likelihood method, making 
use of only the radial velocities of the objects, setting aside 
for the moment the available proper-motion information. 
Specifically, we focus on the difference in the mass esti- 
mate arising from the presence or absence of Leo I. Fig. 

6 shows the likelihood contours in the mass-anisotropy 
(M — (3) plane for the case of a power-law tracer popu- 
lation with 7 = 3.4, where (3 is limited to the range of 
— 1.5 < (3 < 1. The solid and dashed lines denote the 
presence and absence of Leo I, respectively. As is evident, 
the mass estimate sensitively depends on whether or not 
Leo I is bound to the Galaxy, as has been noted in pre- 
vious studies. Inclusion of Leo I yields a likely total mass 
that is an order of magnitude greater than the case with- 
out Leo I. Over the range of (3 we consider, the most likely 
value of M, with Leo I included ,is 21.0 x 10 11 M Q , corre- 
sponding to a scale length a = 185 kpc, whereas excluding 
Leo I yields M = 9.6 x 10 n M Q , and a = 85 kpc. We note 
that the role of Leo I in the Galaxy's mass estimate is also 
understandable from the escape- velocity argument; if only 
the sample radial velocities are taken into account, Leo I 
alone determines the best-fit boundary line Vrf = V csc in 
the Vrf vs. V csc diagram (Fig. 3). 

As is seen in Fig. 6, the high-probability region is bi- 
ased toward the line (3 = —1.5. This bias arises from 
the specific form of the phase-space distribution func- 
tion, F(e,l), given in Eq. (8), where the probability 
P(a, (3\ri,V( r )i,I) is high at large F. We plot F in Fig. 

7 for a set of r and (3 (solid and dotted lines for j3 = — 1 
and 1, respectively). It follows that F at high e is larger 
for smaller (3, whereas F at low e is larger for larger (3. 
The range of e corresponding to these two different cases 
depends on r, as can be deduced from the comparison 
between panel (a) and (b) in Fig. 7. Since our sample 
objects are mainly distributed in the region of higher e 
(solid histograms for the sample with radial velocities), 
the probability is highest at smallest (3. 

Following the above experiments, we drop the lower 
bound of —1.5 for (3, and search for the maximum proba- 
bility at smaller (3. No maximum is found up to (3 = —20, 
although the large discrepancy in M between the cases 
with and without Leo I remains. When we confine our- 
selves to the sample at r > 10 kpc, there exists a maximum 
probability at (3 — —2.75 (with Leo I), with a correspond- 
ing mass 32.0 x 10 11 M Q . For the sample at r > 20 kpc, 
we obtain 11.4 x lO n M at (3 = 0.8. This clearly sug- 
gests that the best-fitting (3, obtained from the analysis 
when only radial velocities are considered, is rather sen- 
sitive to the range of r employed in the sample selection. 
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Fig. 6. (a) Likelihood contours in the plane of the mass, 
M, and velocity anisotropy, (3, obtained from an analysis 
using only radial velocities. The solid and dashed curves 
show the results including Leo I and excluding Leo I, re- 
spectively; the cross and the asterisk show the maxima 
of the probabilities for each case. Contours are plotted at 
heights of 0.32, 0.1, 0.045, and 0.01 of the peak height. 
The spatial distribution of a tracer population is assumed 
to follow a power-law form with 7 = 3.4. (b) Probabilities 
of the mass M at j3 — —1.5, including Leo I (solid line) 
and excluding Leo I (dashed line). 



This in turn affects the number distribution, N(e), which 
is relevant to the likely range of F (Fig. 7). 

With these unavoidable limitations of the present sam- 
ple in mind, Table 2 summarizes the likelihood results for 
the limited range of — 1.5 < (3 < 1, obtained for power-law 
and shadow tracers, using a variety of different priors on 
a and (3. The most likely value of (3 is —1.5 for all cases, 
for the reasons described above. We note that the current 
mass estimate is rather insensitive to the (3 prior. As the 
(3 prior decreases, the estimated mass generally increases, 
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Fig. 7. The distribution function, F, for (3 = — 1 (solid 
lines) and (3 — 1 (dotted lines), at r — 10 kpc (panel a) 
and r — 50 kpc (panel b). Also plotted are the number 
distributions N(e) of the stars when a — 200 kpc, where 
dotted and solid histograms denote the sample with and 
without available proper motions, respectively. The range 
of r for plotting N(e) (r < 20 kpc for panel (a) and 20 < 
r < 80 kpc for panel b) is chosen to approximately match 
that for F. 



and the best-fitting [3 decreases, because the small (3 prior 
is biased toward more tangcntially anisotropic velocity dis- 
tributions than is the large (3 prior. However, since most 
of our sample have high e, the best- fitting (3 remains —1.5 
regardless of whether we adopt the uniform prior or the 
uniform-energy prior for (3. This property makes the mass 
estimate insensitive to the (3 prior. 

Now we apply the Bayesian likelihood method to 
the subsample of objects with both radial velocities and 
proper motions available, and consider the derived space 



motions. In contrast to the above case, where we used ra- 
dial velocities alone, we find that the maximum probabil- 
ity within the range of (3 we consider is now bounded (Fig. 
8a). This may be caused by the characteristic distribution 
of e for the sample with full space motions, as shown in 
Fig. 7 (dotted histogram). This figure shows that there 
exists a larger fraction of low-e stars than are found in 
the sample with radial velocities alone (solid histograms), 
so a larger (3 is preferred to achieve a larger F . The mass 
estimate in this case is quite insensitive to the presence 
or absence of Leo I. Fig. 8b shows the probabilities, as 
a function of M, with a fixed value of (3 = —1.25, for 
the case of a power-law tracer population with 7 = 3.4. 
Solid and dashed lines denote the probabilities with and 
without Leo I, respectively. As is evident, the agreement 
between both probabilities is significantly improved com- 
pared to the case when the radial velocities are considered 
alone (Fig. 6b). When Leo I is included, the most likely 
value of the total mass, M, and the scale length, a, are 
25.0 x lO n M0 and 225 kpc, respectively. Excluding Leo I 
yields M = 18.0 x 10 n M Q and a = 160 kpc. Table 3 
summarizes the various results obtained when the proper 
motions of the objects are considered. This Table illus- 
trates that, for all cases, the mass of the Galaxy obtained 
when including Leo I is in good agreement with that ob- 
tained without Leo I. Also, the mass estimate depends 
only weakly on the index 7, unknown prior probabilities 
for a and (3, as well as on the range of r used in the sample 
selection, resulting in small changes in the mass estimates 
over a range of only a few times 10 11 M Q . 

To estimate the typical errors in this mass determi- 
nation that are associated with the measurement errors 
of the 561 tracers we have analyzed, we have conducted 
Monte Carlo simulations, adopting the assumptions that 
typical errors in the distances and radial velocities are 
10 %, and 10 km s , respectively, and that the proper- 
motion errors are 1 mas yr _1 for globular clusters, 0.3 
mas yr _1 for satellite galaxies, and 5 mas yr -1 for the 
FHB stars. We generated 561 data points (including Leo I) 
drawn from Gaussian distribution functions centered on 
the observational data, and with dispersions set to the 
above typical errors. Given a true mass M, or scale length 
a (where we use M = 2.3 x 10 12 M© with a = 200 kpc), 
and prior probabilities for a and (3 (1/a 2 and the uniform- 
energy prior, respectively), we calculate the most likely 
mass, M', and compare it with an input true mass. Fig. 
9 shows the distribution of the discrepancy between M' 
and M, 100 x (M 1 - M)/M, obtained from 1000 realiza- 
tions. The error distribution in the current mass estimate 
has a mean value shifted downward by 20 %, and a dis- 
persion of half- width 20 %. These values suggest that one 
might adopt an estimate of the systematic error on the 
order of 20 %, and a random error of ±20 %. Exclusion 
of Leo I does not influence the magnitude of these errors. 
It is worth noting that WE99 obtained roughly ~ 100 % 
systematic errors, and ~ 90 % random errors in their mass 
estimate, which was based on about 30 data points. The 
significant improvement of our mass estimate is mainly 
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Fig. 8. (a) Likelihood contours in the plane of the mass, 
M, and velocity anisotropy, /3, obtained from an analysis 
that uses both radial velocities and proper motions. Solid 
and dashed curves show the results including Leo I and ex- 
cluding Leo I, respectively; the cross and the asterisk show 
the maxima of the probabilities for each case. Contours 
are plotted at heights of 0.32, 0.1, 0.045, and 0.01 of the 
peak height. The spatial distribution of a tracer popula- 
tion is assumed to follow a power-law form with 7 = 3.4. 
(b) Probabilities of the mass M at the best-fitting /3 of 
— 1.25, including Leo I (solid line) and excluding Leo I 
(dashed line). 



due to our consideration of a much larger data set that 
includes several hundred FHB stars. 

As shown in Table 3, the most likely estimated total 
mass depends on model assumptions at a level of a few 
times 10 11 Mq. When the model is fixed, the current large 
data set allows us to limit both systematic and random 
errors to a level of about 20 %. If we follow WE99's pro- 
cedure for the adoption of the most likely total mass, i.e., 
if we adopt the mass estimate that provides the smallest 
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Fig. 9. An approximate error distribution of the mass 
estimate caused by the typical measurement errors of the 
data. The abscissa denotes the relative error in mass, 100 x 
(M' — M)/M, where M' is the mass calculated by a Monte 
Carlo method and M is the input true value. See text for 
more details. 

difference between the masses obtained with Leo I and 
without Leo I, we obtain 2.51°^ x 10 12 M© (Leo I in- 
cluded) and 1.8±o 7 x 10 12 M Q (Leo I excluded). On the 
other hand, the mass estimate within the distance of the 
LMC (50 kpc) is quite robust, covering the narrow range 
5.4 to 5.5 xlO n M . 

5. DISCUSSION AND CONCLUDING 
REMARKS 

We have placed new limits on the mass of the Galaxy, 
based on a newly assembled set of halo objects with the 
latest available proper-motion data. First, the comparison 
of their space velocities with the escape velocities at their 
estimated distances allowed us to show that the mass lim- 
its we obtained depend on neither the presence or absence 
of Leo I, nor on the large velocity errors for Draco and 
Pal 3; a possible a lower limit on the total mass of the 
Galaxy is about <~ 2.2 x 1O 12 M . Secondly, a Bayesian 
likelihood approach has been used to derive a total mass 
estimate for the Galaxy that is insensitive to the pres- 
ence or absence of Leo I, at least when proper motions 
are taken into account. Although the best mass estimate 
obtained from this approach depends somewhat on model 
assumptions (prior probabilities for a and [3 and possi- 
bly the shape of F, see below), the resultant systematic 
change of the total mass is confined to within a few times 
10 n MQ. The most likely total mass of the Galaxy we 
derive is 2.5^5'q x 10 12 M@. This is in good agreement 
with the total mass obtained by WE99 (1.9±f f x 10 12 M Q ) 
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and that obtained from other methods (e.g., Peebles 1995, 
2 x 10 12 M Q ). Since the size of our tracer sample is signif- 
icantly larger than used in previous studies, both system- 
atic and random errors are reduced to a great extent. We 
note that consideration of the numerous FHB stars plays 
a vital role in this mass estimate, as demonstrated in Sect. 
3. 

It is also worth noting that, if we fix the mass of 
the Galaxy equal to our most likely mass estimate, there 
is insufficient matter present to gravitationally bind the 
LMC, if we adopt the recent proper- motion measurement 
by Anguita et al. (2000). These authors reported rather 
high proper motions, (/z a cos<5, ^s) — (+1.7 ± 0.2, +2.9 ± 

0. 2), compared to previous measurements, (/i Q cos£, fis) = 
(+1.94 + 0.29,-0.14 + 0.36) (Kroupa & Bastian 1997). 
Thus their results need confirmation from other studies^. 

The current work also implies that the Galactic rota- 
tion curve at outer radii, R > Rq, does not decline out 
to at least R ~ 20 kpc (as long as local disturbances to 
circular motions, such as warping motions and/or non- 
axisymmetric motions, are ignored). As illustrated in Fig. 
2, a declining rotation curve corresponding to a = 20 kpc 
and Vlsr = 211 km s _1 fails to bind many sample objects 
to the Galaxy. The smallest possible value for a to bind 
all objects in the isothermal-like density distribution (Eq. 
4) is a = 195 kpc, yielding Vlsr — 220 km s _1 . 

In a more general context, the detailed shape of the 
rotation curve at and beyond R = Rq reflects the inter- 
play between the disk and halo mass distributions, as this 
region is located near the boundary of both components. 
Thus, determining the rotation curve at Rq 5s R ^ 15 kpc 
will set useful limits on the mass distribution in the inner 
parts of the Galaxy. Indeed, the Japanese project VLBI 
Exploration of Radio Astrometry ( VERA) will be able to 
determine both inner and outer rotation curves from the 
measurement of trigonometric parallaxes and proper mo- 
tions of astronomical maser sources that are widely dis- 
tributed in the Galactic disk (Sasao 1996; Honma et al. 
2000). VERA will reach unprecedented astrometric preci- 
sion, ~ 10/ias, and will yield precise determinations of the 
Galactic constants Rq and Vlsr- We note that whatever 
results are derived for the rotation curve, the total mass 
of the Galaxy ought to be larger than 10 12 M Q , in order 
to bind the more distant stellar objects. 

Our estimate for the mass of the Galaxy inside 50 kpc, 

1. e., within the distance of the LMC, is 5.5^;" x 10 11 M Q 
(Leo I included) and 5.3t^ x 10 11 M (Leo I excluded). 
The error estimates are calculated from the maximum and 
minimum values of the total mass. Thus, about 24% of 
the total mass of the Galaxy resides within r < 50 kpc. 
This implies that the possibility of brown dwarfs as the 
origin of the microlensing events toward the LMC may be 
excluded, because it requires a much smaller mass inside 

2 Since we finished our analysis here, a new study by Pedreros 
et al. (2002) has been published [using the same method as 
Anguita el al. (2000)], that reaches the same conclusion as 
Kroupa & Bastian (1997). 



50 kpc, - 1.3 x 10 n M Q (Honma & Kan-ya 1998). Our 
result is also in good agreement with the recent statistics 
of the microlensing events obtained from analysis of the 
5.7-year baseline of photometry for 11.9 million stars in 
the LMC (Alcock et al. 2000), showing the absence of 
short-duration lensing events by brown dwarfs. However, 
the most recent work has suggested that perhaps one of 
the microlensing events is actually caused by a nearby 
low-mass star in the Galactic disk (Alcock et al. 2001). 
More direct observations for identifying lensing objects 
are required to settle this issue. 

Once the total mass of the Galaxy is fixed, it is possi- 
ble to place a useful constraint on the mass of the Local 
Group. Most of the mass in the Local Group is concen- 
trated in M31 and the Galaxy. The total mass of M31 
can be estimated from the positions and radial veloci- 
ties of its satellite galaxies, globular clusters, and plan- 
etary nebulae (Evans & Wilkinson 2000; Cote et al. 2000; 
Evans et al. 2000). If we take it to be 1.2±J;| x 10 12 M Q 
(Evans & Wilkinson 2000), the mass of the Local Group 
is ~ 3.7 x 1O 12 M . This is in good agreement with the 
estimate by Schmoldt & Saha (1998), (4 - 8) x 10 12 M Q , 
based on modified variational principles. 

To set tighter limits on the total mass of the Galaxy 
we require more accurate proper-motion measurements 
for a greater number of objects at large Galactocentric 
distances. The high-velocity FHB stars in our sample 
(with apparent magnitudes V < 16) that are responsi- 
ble for setting the minimum mass of the Galaxy have 
proper- motion errors of ~ 5 mas yr _1 , whereas Draco 
and Pal 3 have much larger relative errors, comparable 
to their proper motions themselves (see Table 2). Indeed, 
both the Space Interferometry Mission (SIM: Unwin et 
al. 1997) and the Global Astrometry Interferometer for 
Astrophysics (GAIA: Lindegren & Perryman 1996) will 
be able to provide more accurate proper motions for such 
high-velocity objects, as well as for numerous other dis- 
tant tracers of the Galaxy's mass, up to a precision of 
a few /Ltas for targets with V < 15. This corresponds 
to an error of ^ 10 km s -1 in the tangential velocity 
components for many distant objects, i.e., comparable to 
the error of their (presently determined) radial veloci- 
ties. Furthermore, roughly half of our sample objects lack 
proper- motion measurements altogether. To a great ex- 
tent, the lack of proper-motion measurements (at least for 
southern sources) will be removed with the completion of 
the recently re-started Southern Proper Motion survey of 
van Altena and colleagues, as well as other efforts to sub- 
stantially increase the numbers of stars with reasonably 
well-measured proper motions (e.g., UCAC1: Zacharias et 
al. 2000; UCAC2: Zacharias et al. 2001). 

Further assembly of radial velocities for FHB stars, es- 
pecially those at large r (beyond distances where accurate 
ground-based proper motions can be obtained), is also of 
great importance for a number of reasons. First, as Fig. 
3 demonstrates, large Galactocentric regions are charac- 
terized by small escape velocities. The current sample of 
FHB stars (because of their locations near the Sun) ex- 
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plore distances where the corresponding escape velocities 
are in the range of 500 ^ U csc ^ 600 km s -1 . More re- 
mote FHB stars, with distances in the range 10 r & 50 
kpc, will offer a further constraint on the total mass of the 
Galaxy by covering the range 400 ;$ V csc ^ 500 km s _1 . 
Secondly, the assembly of samples of more distant FHB 
stars will enable exploration of the suggested change in 
velocity anisotropy from the inner to the outer halo (e.g., 
Sommer-Larsen et al. 1997), and better constrain its de- 
pendence on Galactocentric distance. 

In exploring the Bayesian approach for mass estimates 
of the Galaxy, we have adopted a specific form of the 
phase-space distribution function F (Eq. 8) to facilitate 
comparison with previous studies. This procedure implic- 
itly assumes that the velocity-anisotropy parameter, /3, is 
constant everywhere in the Galactic volume. However, as 
noted by Sommer-Larsen et al. (1997), there is an indi- 
cation that the velocity anisotropy of the halo may be 
mostly radial at R 20 kpc and tangential at R 20 
kpc. If so, many of distant FHB stars, especially those at 
R > 20 kpc, play a crucial role in the determination of 
the global distribution of velocity anisotropy. Searches for 
a more realistic form of the phase-space distribution func- 
tion, combined with a more elaborate likelihood method, 
are both worthy pursuits. Also, instead of exploring such a 
specific but realistic form of distribution function, a non- 
parametric method as proposed by Merritt & Trcmblay 
(1993) will be more useful if a much larger data set is 
available. 

Moreover, the implicit assumption behind the cur- 
rent method, that the sample stars have random distri- 
butions in location and in their space motions, may not 
be well satisfied if the halo is largely dominated by co- 
herent structures such as tidal streams (e.g., Ibata et al. 
2001). Alternative approaches to obtaining mass limits us- 
ing tidal streams (Johnston et al. 1999) are worth consid- 
ering in such a case. 

Fortunately, prospects are excellent for obtaining a 
rapid increase in the observational database of FHB stars 
with the required data. There already exists a substan- 
tial body of additional spectroscopy for FHB/ A stars ob- 
served during the course of the HK survey of Beers and col- 
leagues and the Hamburg/ESO Stellar survey (Christlieb 
et al. 2001), many of which also have available proper 
motions, or will soon, from completion of the SPM sur- 
vey and/or other ground-based efforts. However, as was 
noted by Wilhelm et al. (1999a) (foreshadowed by Norris 
& Hawkins 1991; Rodgers & Roberts 1993, and references 
therein; Kinman et al. 1994; Preston et al. 1994), a sub- 
stantial fraction (perhaps as high as 50%) of high-latitude 
A-type stars are not FHB, but rather some (as yet unde- 
termined) mixture of binaries and high-gravity stars (see 
Preston & Sneden 2000). For some applications, such as 
estimates of the mass of the Galaxy that rely on space mo- 
tions of tracers (and in turn on reasonably precise distance 
estimates of individual objects), confident separation of 
bona-fide members of the FHB population from possible 



"contaminants" is crucial []. In the past, this has required 
that one obtain either Stromgren photometry and/or spec- 
trophotometry (e.g., Kinman et al. 1994), broad-band 
UBV photometry in combination with medium-resolution 
spectroscopy (e.g., Wilhelm et al. 1999a), or reasonably 
high S/N, high-resolution spectroscopy (e.g., Preston & 
Sneden 2000). All such endeavors are rather time inten- 
sive. However, Christlieb et al. (2002, priv. comm.) have 
been exploring means by which adequate separation of 
FHB stars from higher-gravity A-type stars might be ac- 
complished directly from objective-prism spectra, such as 
those in the Hamburg/ESO stellar survey. Such methods, 
which look promising, would be most helpful in future in- 
vestigations of this sort. Wide-field stellar surveys, such as 
those presently being carried out with the 6dF facility at 
the UK Schmidt Telescope, are capable of providing large 
numbers of radial velocities for FHB/A candidates, and 
are expected to contribute 5,000-10,000 suitable data over 
the course of the next few years. 
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Table 2. Likelihood Results for Only the Radial Velocities 



a s or 7 a prior (3 prior Leo I best (3 best a (kpc) best M a M(< 50kpc) a M(< 100kpc) a 



Power-law Tracers 



7 = 3.4 


1/a 2 


Energy 


Yes 


-1 160 


18.0 


5.4 


9.6 








No 


-1 70 


7.9 


4.6 


6.5 


7 = 3.4 


1/a 


Energy 


Yes 


-1 175 


20.0 


5.4 


9.8 








No 


-1 75 


8.5 


4.7 


6.8 


7 = 3.4 


1/a 2 


Uniform 


Yes 


-1 160 


18.0 


5.4 


9.6 








No 


-1 70 


7.9 


4.6 


6.5 


7 = 4.0 


1/a 2 


Energy 


Yes 


-1 170 


19.0 


5.4 


9.7 








No 


-1 80 


9.0 


4.8 


7.1 


Shadow Tracers 


a s = 100 


1/a 2 


Energy 


Yes 


-1 180 


20.0 


5.4 


9.8 








No 


-1 65 


7.4 


4.5 


6.2 


a s = 100 


1/a 


Energy 


Yes 


-1 205 


23.0 


5.5 


10.0 








No 


-1 70 


7.9 


4.6 


6.5 


a s = 100 


1/a 2 


Uniform 


Yes 


-1 180 


20.0 


5.4 


9.8 








No 


-1 65 


7.4 


4.5 


6.2 


&s — ^halo 


1/a 2 


Energy 


Yes 


-1 170 


19.0 


5.4 


9.7 








No 


-1 70 


7.9 


4.6 


6.5 


a All masses 


are in 


units of 10 11 


M Q . 














Table 3. 


Likelihood Results for the Full Space Velocities 




a s or 7 


a prior fi prior 


Leo I 


best (3 best a (kpc) 


best M a 


M(< 50kpc) a 


M(< 100kpc) a 



Power-law Tracers 



7 = 3.4 


1/a 2 


Energy 


Yes 


-1 200 


23.0 


5.5 


10.0 








No 


-1 150 


17.0 


5.3 


9.4 


7 = 3.4 


1/a 


Energy 


Yes 


-1 230 


26.0 


5.5 


10.0 








No 


-1 160 


18.0 


5.4 


9.6 


7 = 3.4 


1/a 2 


Uniform 


Yes 


-1 200 


23.0 


5.5 


10.0 








No 


-1 150 


17.0 


5.3 


9.4 


7 = 4.0 


1/a 2 


Energy 


Yes 


-1 225 


25.0 


5.5 


10.0 








No 


-1 165 


19.0 


5.4 


9.6 


Shadow Tracers 


a s = 100 


1/a 2 


Energy 


Yes 


-1 275 


31.0 


5.5 


11.0 








No 


-1 190 


21.0 


5.4 


10.0 


a s = 100 


1/a 


Energy 


Yes 


-1 330 


37.0 


5.6 


11.0 








No 


-1 220 


25.0 


5.5 


10.0 


a s = 100 


1/a 2 


Uniform 


Yes 


-1 275 


31.0 


5.5 


11.0 








No 


-1 190 


21.0 


5.4 


10.0 


fl s — Q-halo 


1/a 2 


Energy 


Yes 


-1 275 


31.0 


5.5 


11.0 








No 


-1 190 


21.0 


5.4 


10.0 



All masses are in units of 10 11 M Q . 



TABLE 1 . Distance Estimates and Proper Motions of the FHB stars 



NAME* 




RA 


DEC 


[Fe/Hl 


r b 


RV C 


fta * 




V 


Source d 




(2000.0) 








(km/s) 


(mas/ vr ) 


(mas / vr) 


(km/s) 






228760029 





:01:57.6 


-36:40:46 


-2.8 


8.8 


48 


-0.8 ± 3.4 


-5.6 ± 2.8 


148 ± 90 




S 


295170031 





:02:09.8 


-14:08:55 


-2.7 


12.4 


-314 


5 ± 5 


-10 ± 5 


394 ± 211 




N 


295170044 





;03:08.4 


-14:24:25 


-2.6 


10.4 


60 


-9 ± 5 


-1 ± 5 


477 ± 163 




N 


228760030 





:05:36.5 


-36:41:28 





7.9 


50 


2.0 ± 2.4 


1.7 ± 2.1 


244 ± 28 




A, S 


295030008 





:06:02.2 


-24:37:09 


-2 


9.3 


45 


1.4 ± 3.4 


-3.0 ± 4.9 


156 ± 120 




S 


228760034 





:06:20.7 


-35:17:14 


-1.2 


8.2 


-94 


3.3 ± 10.0 


-6.4 ± 8.9 


140 ± 127 




S 


228760031 





:08:25.6 


-36:09:15 


-2.2 


8.7 


-54 


-2.2 ± 2.8 


-8.7 ± 2.8 


193 ± 62 




S 


228760038 





:09:44.7 


-34:39:14 


-1.9 


8.8 


-153 


4.4 ± 2.2 


-5.4 ± 2.3 


177 ± 31 




S 


295030024 





:10:08.5 


-25:33:40 


-1.4 


9.7 


143 


13.2 ± 3.6 


-3.7 ± 2.9 


347 ± 87 




s 


295030029 





;11:19.1 


-26:26:39 


-2.4 


8.6 


-49 


12.7 ± 2.0 


-16.6 ± 1.8 


197 ± 44 




A, S 


295270016 





:27:53.9 


-18:57:44 


-3 


9.5 


-43 


11 ± 5 


-3 ± 5 


193 ± 125 




N 


294970009 





:29:01.7 


-23:40:01 


-2.1 


9.3 


-131 


15.5 ± 3.9 


-2.4 ± 5.1 


308 ± 100 




S 


295270022 





:32:17.1 


-19:26:00 


-1.9 


8.2 


-98 


25.0 ± 2.0 


-8.1 ± 2.0 


168 ± 18 


A, 


N, T 


295270026 





:33:05.0 


-21:11:13 


-1.8 


9.4 


-146 


6 ± 5 


-3 ± 5 


176 ± 76 




N 


221700002 





:34:37.9 


-10:28:50 


-2.5 


10.0 


-270 


-4 ± 5 


-10 ± 5 


321 ± 102 




N 


295270035 





:35:23.2 


-21:00:47 


-2.2 


8.7 


-22 


13.9 ± 3.8 


-11.9 ± 3.8 


97 ± 63 




A, N 


228820023 





;36:04.3 


-30:16:11 


-1.8 


10.6 


-132 


-2.3 ± 2.7 


-9.1 ± 2.8 


291 ± 83 




S 


295270031 





:36:39.2 


-22:25:44 


-2.1 


9.1 


-19 


3.7 ± 2.2 


-13.1 ± 2.9 


92 ± 51 


A, 


N, S 


295270039 





:36:53.8 


-19:56:51 


-1.1 


10.0 


-26 


2 ± 5 


2 ± 5 


268 ± 138 




N 


294970031 





:38:41.9 


-24:26:56 


-2 


9.3 


-24 


5.8 ± 2.5 


-2.0 ± 2.5 


148 ± 57 




S 


221790011 





;39:36.1 


-04:35:40 


-3 


9.1 


-110 


3.0 ± 4.4 


-10.0 ± 4.4 


95 ± 60 




A 


221700009 





:40:52.9 


-10:19:31 


-2.2 


9.9 


-125 


9 ± 5 


-8 ± 5 


105 ± 92 




X 


294970033 





:41:08.5 


-25:56:31 


-1.9 


9.2 


73 


1.6 ± 4.0 


-9.4 ± 2.7 


124 ± 65 




s 


221700013 





;41:21.1 


-08:23:19 


-2 


9.1 


-7 


16.8 ± 3.3 


-30.9 ± 3.3 


50 ± 56 




A, X 


221700015 





:43:01.6 


-08:15:10 





9.6 


-22 


-11 ± 5 


-13 ± 5 


362 ± 100 




X 


294970038 





:43:47.1 


-26:43:50 


-3 


13.1 


-42 


0.0 ± 4.8 


-6.5 ± 4.1 


192 ± 222 




S 


221700024 





;43:47.3 


-11:22:15 


-2.2 


12.0 


81 


2 ± 5 


-6 ± 5 


150 ± 105 




X 


295270061 





:43:48.6 


-20:45:41 


-3 


9.5 


-67 


7.0 ± 4.1 


-3.9 ± 4.1 


118 ± 81 




A, X 


221830011 





:52:56.1 


-02:52:42 


-1.7 


8.9 


-126 


13.9 ± 2.1 


-31.5 ± 2.1 


187 ± 25 




A, T 


295090031 





:52:57.0 


-29:54:38 


-1.6 


9.3 


131 


13.5 ± 2.4 


-11.0 ± 1.9 


219 ± 44 




S 


295090039 





:54:41.4 


-28:13:54 


-1.7 


9.0 


54 


15.8 ± 2.1 


-7.4 ± 2.2 


186 ± 42 




s 


221830014* 


1 


:00:15.3 


-02:17:29 


-3 


10.9 


85 


2 ± 5 


16 ± 5 


635 ± 127 




X 


221830028 


1 


:02:07.6 


-07:02:34 


-2.2 


10.1 


-50 


2 ± 5 


12 ± 5 


441 ± 101 




X 


221830024 


1 


:02:53.2 


-04:44:54 


-2.2 


8.7 


-100 


11.2 ± 2.1 


-15.5 ± 2.1 


70 ± 16 




A, T 


295140008 


1 


;05:42.6 


-23:58:40 


-1.3 


9.0 


113 


4.8 ± 3.0 


-8.7 ± 2.3 


129 ± 25 




A, S 


295140006 


1 


:06:49.4 


-25:01:48 


-1 


9.8 


57 


5.5 ± 5.7 


-10.2 ± 5.7 


77 ± 104 




S 


221660032 


1 


:07:33.8 


-12:11:01 


-2.6 


9.5 


-128 


15.0 ± 3.1 


-7.3 ± 3.1 


166 ± 41 




A, X 


221660034* 


1: 


;08:03.3 


-12:11:37 


-2.2 


11.2 


128 


4 ± 5 


13 ± 5 


590 ± 149 




X 


295140013 


1 


:10:31.3 


-25:40:47 


-2.7 


10.3 


-50 


0.8 ± 3.8 


-12.9 ± 3.7 


226 ± 105 




S 


295180028 


1 


:16:44.5 


-31:04:58 


-2.5 


11.7 


14 


0.5 ± 2.5 


-4.2 ± 2.6 


127 ± 100 




S 


295180035 


1 


:16:58.2 


-27:45:57 


-2.4 


10.7 


45 


-3.4 ± 4.3 


-9.8 ± 4.3 


280 ± 133 




S 


295140038 


1 


;22:53.5 


-26:17:35 


-1.6 


9.8 


54 


5.9 ± 4.6 


-14.3 ± 3.1 


154 ± 87 




S 


221740034 


1 


:25:24.5 


-09:36:19 


-2.5 


9.8 


37 


1.6 ± 3.4 


-10.2 ± 3.4 


136 ± 51 


A, 


X, T 


221740042 


1 


:30:19.3 


-09:44:57 


-2.5 


9.0 


7 


4.3 ± 1.8 


-18.6 ± 1.8 


109 ± 19 


A, 


X, T 



Table 1. (continued) 



IN t\ IV 1 JZj 




RA 


DEC 




r b 


RV C 


Met* 






o d 

oourcG 




(2000.0) 






(KpC) 


{KILL/ o J 


(mas/yr) 


(mas/ yr) 






221800003 


1 


:30:46.2 


-10:25:54 


-2.7 


10.6 


79 


9 + 5 


11 + 5 


442 ±114 


X 


295040004 


1 


:31:02.1 


-36:38:34 


-2.3 


10.1 


45 


8.6 ± 2.4 


-5.6 ± 2.6 


100 ± 68 


s 


221800002 


1 


:31:42.6 


-10:05:30 


-1.1 


9.3 


-82 


10 + 5 


3 + 5 


232 ± 60 


X 


221800006 


1 


:35:19.5 


-10:33:20 


-2.2 


11.0 


47 


5 + 5 


6 + 5 


350 ± 128 


N 


295040028 


1 


:36:12.2 


-34:06:08 


-2.3 


10.0 


-56 


6.7 ± 3.8 


-8.4 ± 3.4 


120 ± 32 


S 


221800017 


1 


:36:49.9 


-12:00:53 


-1.5 


9.6 


34 


11.1 ± 3.1 


-5.2 ± 3.1 


111 ± 50 


A, N 


295040035 


1 


:40:51.5 


-33:25:36 


-2.2 


11.1 


43 


10.6 ± 4.8 


-12.0 ± 4.8 


310 ± 156 


S 


295040045 


1 


:47:30.5 


-34:07:23 


-2 


9.8 


189 


9.0 ± 3.8 


-6.0 ± 2.4 


144 ± 49 


S 


221710022 


2 


:03:11.7 


-08:13:10 


-2.2 


9.7 


-201 


-2.4+ 2.8 


-22.7 ± 2.8 


310 ± 38 


A, N 


221750001 


2 


:13:30.0 


-11:37:36 


-1.4 


9.7 


67 


5.1 ± 2.9 


4.2 ± 2.9 


256 ± 41 


A, N 


221750003 


2 


;15:32.1 


-10:40:28 


-1.8 


10.7 


16 


14 + 5 


6 + 5 


335 ± 91 


X 


221890005 


2 


:32:29.2 


-14:31:48 


-1.6 


11.4 


-120 


10 + 5 


— 5 ± 5 


177 ± 60 


X 


221810032 


3 


;01:34.1 


-08:56:03 


-1.2 


11.5 


-126 


5 + 5 


2 + 5 


273 ± 96 


X 


310640037 


3 


;09:30.0 


-67:37:08 


-1.5 


7.9 


277 


67.1 ± 2.1 


-41.8 ± 2.1 


496 ± 22 


A, T 


221670008 


3 


:16:28.3 


-07:06:46 





15.2 


14 


1 ± 5 


+ 5 


204 ± 209 


X 


221670017 


3 


:18:45.7 


-03:38:50 


-2 


11.2 


198 


14.0 ± 3.9 


-21.0 ± 3.9 


322 ± 74 


A 


221850020 


3 


;27:47.2 


-14:06:46 


-1.7 


9.7 


24 


19.7 ± 3.1 


-22.8 ± 3.1 


144 ± 37 


A, X 


310750042 


3 


:30:46.9 


-66:38:17 


-2.1 


8.2 


145 


13.5 ± 4.4 


1.0 ± 4.4 


86 ± 67 


A 


221760020 


3 


:45:49.2 


-10:23:13 


-1.5 


10.9 


-6 


16 + 5 


-14 ± 5 


181 ± 81 


X 


221690025 


4 


;16:05.6 


-14:34:02 


-1.6 


12.1 


10 


17+5 


2 + 5 


336 ± 93 


X 


310720061 


5 


:27:10.8 


-59:05:17 


-1.5 


8.5 


64 


-0.5 ± 6.6 


-1.0 ± 6.6 


188 ± 57 


A 


310720068 


5 


:29:20.6 


-61:16:28 


-2.4 


8.2 


495 


41.8 ± 3.4 


-1.8 ± 3.4 


369 ± 20 


A, T 


156210043 


10 


:09:27.2 


+24:50:05 


-3 


12.5 


67 


4.5 ± 4.5 


-2.0 ± 4.5 


245 ± 113 


A 


156210039 


10 


:13:50.7 


+ 25:18:24 





10.3 


131 


-4.5 ± 5.0 


-8.0 ± 5.0 


113 ± 71 


A 


156210070 


10 


:21:56.1 


+ 27:11:19 


-1.8 


9.6 


35 


17.0 ± 1.8 


-6.5 ± 1.8 


304 ± 20 


A, T 


156210015 


10 


:26:19.1 


+23:30:36 


-2.2 


9.3 


181 


-23.7+ 1.4 


-11.4 ± 1.4 


246 ± 17 


A, T 


156210009 


10 


:28:01.0 


+ 27:32:38 


-1.4 


11.0 


27 


-17.7 ± 5.0 


-2.0 ± 5.0 


373 ± 91 


A 


156210010 


10 


;28:06.1 


+26:32:51 





10.5 


46 


0.0 ± 7.2 


-28.0 ± 7.2 


324 ± 138 


A 


156250026 


11 


:51:50.0 


+26:21:56 


-1.1 


8.3 


-3 


7.8 ± 1.6 


-17.8 ± 1.6 


189 ± 10 


A, N, T 


160260011 


12 


:16:50.0 


+ 28:56:03 


-1.7 


9.3 


19 


-9.7 ± 3.3 


-8.5 ± 3.3 


70 ± 58 


A, X 


160260028 


12 


;23:02.8 


+27:27:15 


-3 


9.6 


-81 


5.2 ± 3.5 


8.9 ± 3.5 


467 ± 77 


A, N 


160270049 


13 


:12:26.9 


+30:21:16 


-1.9 


8.7 


57 


-14.4 ± 3.6 


-10.3 ± 3.6 


143 ± 49 


A, N 


228770008 


13 


:12:49.9 


-10:32:31 


-2.2 


7.8 


90 


+ 5 


-5 + 5 


135 ± 140 


X 


160270051* 


13 


:13:01.9 


+31:01:28 


-2.2 


11.8 


39 


-6 + 5 


8 + 5 


557 ± 211 


X 


228770012 


13 


;13:32.9 


-09:35:18 


-2.6 


7.4 


312 


4.0 ± 3.7 


0.6 ± 3.7 


355 ± 57 


A, N 


228770005 


13 


:15:28.5 


-11:25:31 


-3 


8.4 


31 


2 + 5 


1 ± 5 


308 ± 164 


X 


228770031 


13 


:17:08.5 


-09:40:25 


-2.7 


7.6 


65 


-6.8 ± 5.0 


2.3 ± 5.0 


225 ± 104 


A, N 


228770036 


13 


:17:39.4 


-11:18:33 


-2.3 


7.6 


202 


-8.7 ± 5.0 


1.2 ± 5.0 


239 ± 99 


A, N 


228770038 


13 


:18:01.9 


-11:58:46 


-2.2 


7.2 


57 


-13.4 ± 5.0 


4.3 ± 5.0 


270 ± 82 


A, X 


228770027 


13 


:18:10.7 


-08:51:05 


-1.9 


8.6 


-11 


-7.4 ± 5.0 


-10.6 ± 5.0 


292 ± 147 


A, X 


228770030 


13 


:18:26.0 


-09:07:12 


-3 


9.5 


136 


-7.8 ± 5.0 


-1.7 ± 5.0 


227 ± 204 


A, N 


228770020 


13 


:18:32.8 


-07:49:17 


-2.2 


7.5 


-18 


-10 ± 5 


4 + 5 


277 ± 99 


X 


228770026 


13 


:20:55.2 


-08:38:48 


-2 


7.3 


161 


-10.3 ± 3.3 


-23.9 ± 3.3 


298 ± 55 


A, N 


228770045 


13 


:21:35.2 


-11:29:15 


-2.1 


7.2 


42 


4.9 ± 5.0 


-16.8 ± 5.0 


272 ± 75 


A, N 



Table 1. (continued) 



NAME" RA DEC [Fe/H] r b RV C fi a . n 6 v Source d 

(2000.0) (kpc) (km/s) (mas/yr) (mas/yr) (km/s) 



228770046 


13:22:24.0 


-11:17:27 


-2.3 


7.3 


60 


-5.6 ± 5.0 


-5.8 ± 5.0 


37 ± 96 


A, N 


228770049 


13:22:52.2 


-11:58:37 


-2 


7.2 


-7 


-10.4 ± 5.0 


-11.4 ± 5.0 


146 ± 78 


A, N 


228890013 


13:36:32.1 


-08:15:15 


-1.9 


7.1 


32 


-5 + 5 


-3 + 5 


129 ± 90 


N 


228890023 


13:38:14.8 


-12:09:25 


-1.9 


7.2 


-75 


-11 ± 5 


-9 + 5 


256 ± 126 




228890037 


13:40:14.5 


-09:14:04 


-2.5 


7.0 


174 


-6 + 5 


-14 ± 5 


173 ± 97 


X 


228890040 


13:42:53.8 


-07:49:02 


-3 


8.9 


-80 


+ 5 


-4 + 5 


220 ± 182 


N 


228890060 


13:49:13.8 


-10:22:21 


-1.1 


6.9 


-107 


-5.9 ± 3.5 


-6.5 ± 3.5 


206 ± 42 


A, N 


228890058 


13:51:13.4 


-11:03:31 





6.8 


65 


6.2 ± 3.2 


-24.3 ± 3.2 


351 ± 46 


A, N 


228890057 


13:52:21.3 


-11:07:54 


-0.6 


6.7 


8 


-11 ± 5 


-2 + 5 


158 ± 101 


N 


156230001 


14:08:30.2 


+ 23:20:51 


-2.8 


8.1 


-159 


11.0 ± 4.7 


-12.0 ± 4.7 


431 ± 89 


A 


228830007 


14:15:34.5 


+ 10:21:43 


-2.6 


7.3 


51 


-20.5 ± 1.8 


-39.0 ± 1.8 


304 ± 35 


A, X, T 


228830028 


14:20:47.1 


+09:15:11 


-3 


7.1 


-1 


-29.6 ± 7.3 


-14.0 ± 7.3 


462 ± 145 


A 


228740025* 


14:30:24.7 


-24:00:11 


-2.1 


6.0 


-27 


-41.1 ± 6.4 


9.0 ± 6.4 


585 ± 85 


A 


228710008 


14:32:03.4 


-20:34:17 


-2 


5.6 


83 


-1 + 5 


-8 + 5 


147 ± 132 


X 


228740021 


14:33:17.2 


-24:32:46 


-2.3 


5.7 


91 


-13.7 ± 6.0 


-20.0 ± 6.0 


278 ± 98 


A 


228710003 


14:33:57.1 


-21:12:58 


-1.5 


5.5 


78 


-16 ± 5 


-3 + 5 


281 ± 113 


X 


228710013 


14:34:14.4 


-19:37:26 


-3 


5.8 


-14 


-8 + 5 


-9+ 5 


204 ± 157 


N 


228710005 


14:34:30.7 


-20:50:56 


-2.2 


6.1 


105 


-11 ± 5 


+ 5 


318 ± 168 


X 


228740036 


14:34:57.5 


-24:13:04 


-1.9 


6.0 


40 


-25.5 ± 6.1 


-13.0 ± 6.1 


260 ± 79 


A 


228710009 


14:35:30.9 


-20:27:53 





7.0 


-179 


-15.7 ± 2.1 


7.7 ± 2.1 


333 ± 15 


A, X, T 


228740053 


14:35:37.4 


-26:31:15 


-0.6 


5.9 


-94 


-25.1 ± 5.8 


0.0 ± 5.8 


339 ± 76 


A 


228710045 


14:38:12.2 


-21:45:48 


-2.3 


5.4 


-5 


-11 ± 5 


-0.2 ± 5 


226 ± 112 


X 


228710031 


14:38:47.6 


-18:54:01 


-3 


5.5 


210 


-12 ± 5 


-7+5 


241 ± 114 


X 


228710034 


14:39:20.1 


-19:16:13 


-1.5 


6.1 


84 


-14 ± 5 


-13 ± 5 


530 ± 198 


X 


228710077 


14:41:17.3 


-18:03:25 


-1.7 


6.1 


-10 


-10 ± 5 


3 + 5 


354 ± 148 


X 


228710062 


14:42:05.6 


-20:38:44 


-2 


5.5 


35 


-5 + 5 


-5 + 5 


35 ± 116 


X 


228710064 


14:42:22.4 


-20:31:01 


-2 


5.4 


192 


4 + 5 


-12 ± 5 


352 ± 121 


X 


228710052 


14:42:53.2 


-22:24:30 


-2.3 


5.5 


79 


-8 + 5 


-1 + 5 


130 ± 88 


X 


228710063 


14:43:00.6 


-20:42:53 


-3 


5.6 


159 


-20 ± 5 


-6 + 5 


259 ± 83 


X 


228710057 


14:43:10.0 


-21:11:28 


-1.9 


5.8 


82 


-16 ± 5 


-8 + 5 


124 ± 71 


X 


228710085 


14:43:51.6 


-18:43:19 


-2.2 


5.4 


209 


-16 ± 5 


-2 + 5 


302 ± 104 


X 


228710088 


14:44:26.0 


-19:20:25 


-1.8 


5.4 


115 


-9 + 5 


3 + 5 


245 ± 105 


X 


228710087 


14:45:16.9 


-18:58:35 


-2.4 


5.4 


-107 


-9 + 5 


11 + 5 


464 ± 103 


X 


228710092 


14:46:08.5 


-20:30:13 


-2.4 


5.3 


57 


-8 + 5 


-5 + 5 


45 ± 105 


X 


228710096 


14:48:02.9 


-21:32:53 


-2 


5.6 


-40 


-17 ± 5 


-13 ± 5 


208 ± 75 


X 


228710115 


14:48:35.2 


-17:37:22 


-2.7 


5.5 


17 


+ 5 


6 + 5 


353 ± 99 


X 


228710109 


14:48:42.3 


-18:47:59 


-2.6 


5.3 


-54 


+ 5 


3 + 5 


302 ± 111 


X 


228710099 


14:50:24.1 


-21:43:10 


-2.5 


5.3 


162 


-3 + 5 


-9 + 5 


190 ± 120 


X 


228710103 


14:50:52.9 


-20:53:17 


-1.2 


5.4 


75 


1 ± 5 


-7 + 5 


209 ± 156 


X 


228710113 


14:51:26.5 


-19:31:48 


-2.8 


5.3 


97 


+ 5 


-11 ± 5 


270 ± 135 


X 


228900015 


15:16:10.2 


+02:09:16 


-3 


6.5 


129 


-25.0 ± 4.3 


-40.0 ± 4.3 


380 ± 60 


A 


228900042 


15:18:53.8 


+00:45:52 


-2.2 


6.4 


-78 


-21.0 ± 5.0 


-13.0 ± 5.0 


130 ± 55 


A 


228840015 


15:31:46.9 


-09:28:32 


-2.1 


5.9 


223 


-4.8 ± 3.3 


-9.4+ 3.3 


241 ± 32 


A, X 


228840021 


15:34:44.9 


-08:16:39 


-1.6 


5.0 


-59 


-13 ± 5 


6 + 5 


358 ± 92 


X 



Table 1. (continued) 



NAME" 


RA 


DEC 


[Fe/H] 


r b 


RV C 






Us 


V 


Source d 




(2000.0) 






(kpc) 


(km/s) 


(mas 


TV) 


(mas/yr) 


(km/s) 




228840006 


15:35:13.9 


-11:40:54 





4.8 


34 




7 + 5 


-6 + 5 


48 ± 106 


N 


228840036 


15:37:38.5 


-11:30:38 


-2.2 


5.6 


69 


-13.7+ 3.5 


-13.1 ± 3.5 


88 ± 49 


A, N 


228840047 


15:41:27.7 


-11:28:06 


-2.1 


5.7 


-90 


-14.1 


± 3.1 


1.6 ± 3.1 


225 ± 36 


A, N 


228720040 


16:21:43.7 


-03:08:00 


-1.2 


5.0 


-65 




2 + 5 


3 + 5 


295 ± 94 


N 


228720041 


16:22:40.4 


-03:24:38 


-2 


6.6 


-204 


-48.9 


± 3.9 


-31.0 ± 3.9 


317 ± 33 


A 


228720071 


16:25:40.3 


-02:57:18 


-3 


4.3 


74 


3.0 


± 8.5 


-3.0 ± 8.5 


278 ± 214 


A 


228720067 


16:25:51.5 


-03:38:59 


-2.6 


4.3 


-16 


1.0 


± 5.8 


-8.0 ± 5.8 


180 ± 134 


A 


228780105 


16:50:04.4 


+08:11:15 


-2.3 


5.7 


-179 


-4.0 


± 5.2 


14.0 ± 5.2 


416 ± 72 


A 


229590022 


18:45:25.7 


-65:57:31 


-1.7 


6.1 


304 


-22.8 


± 4.3 


14.0 ± 4.3 


481 ± 45 


A 


229590189 


19:10:50.0 


-66:33:16 


-1.6 


6.0 


-129 


-0.8 


± 4.8 


-9.0 ± 4.8 


248 ± 54 


A 


229390058 


19:22:28.0 


-28:09:09 


-2 


3.2 


17 


-2.5 


± 3.3 


3.2 ± 3.4 


322 ± 83 


S 


229390167 


19:28:43.8 


-28:50:26 


-2 


5.0 


-32 


-2.7 


± 8.7 


3.3 ± 3.6 


277 ± 71 


S 


228960041 


19:29:22.3 


-52:54:06 


-3 


5.2 


-122 


-0.6 


± 5.8 


-11.0 ± 5.8 


180 ± 76 


A 


229390164 


19:32:31.4 


-28:41:21 


-1.7 


5.3 


-264 


13.7 


± 2.8 


-13.8 ± 2.8 


336 ± 25 


A, S 


228960086 


19:36:45.3 


-57:18:59 





5.5 


-169 


-11.9 


± 5.7 


-9.0 ± 5.7 


284 ± 83 


A 


229640125 


19:57:17.9 


-39:01:18 


-2.1 


4.2 


26 


1.3 


± 2.9 


3.9 ± 2.9 


338 ± 72 


S 


229640219 


20:06:29.4 


-39:03:41 


-0.3 


4.6 


19 


-5.6 


± 2.9 


-1.1 ± 2.9 


230 ± 61 


S 


229500023 


20:16:20.7 


-14:33:38 


-1.3 


7.0 


-13 


19.3 


± 2.0 


-14.4 ± 2.0 


208 ± 12 


T 


229500008 


20:18:25.5 


-16:01:17 


-1.2 


6.0 


-6 


-2.9 


± 5.4 


-3.0 ± 5.4 


196 ± 59 


A 


228850092 


20:26:02.3 


-40:44:34 


-1.8 


4.6 


36 


-1.1 


± 2.3 


-7.0 ± 2.3 


60 ± 51 


S 


228850125 


20:26:53.6 


-37:57:37 


-1.1 


4.8 


-175 


1.6 


± 4.1 


-8.5 ± 4.9 


171 ± 59 


S 


229550104 


20:31:35.1 


-25:51:22 


-1.6 


5.2 


-122 


5.0 


± 2.7 


-4.3 ± 2.6 


230 ± 89 


s 


229550103 


20:31:54.4 


-26:06:51 


-1.6 


7.7 


-180 


2.5 


± 3.0 


-6.9 ± 3.0 


268 ± 139 


s 


229550099 


20:32:47.2 


-26:33:23 


-1.1 


6.5 


-101 


9.4 


± 2.4 


-25.9 ± 2.6 


109 ± 18 


s 


228850179 


20:33:41.1 


-38:42:26 


-1.2 


4.7 


-28 


-0.3 


± 2.2 


0.4 ± 2.2 


246 ± 65 


s 


229550147 


20:38:35.6 


-25:48:06 


-1.1 


5.2 


10 


-6.8 


± 2.0 


-5.8 ± 2.0 


188 ± 38 


s 


228800061 


20:41:28.9 


-21:47:13 


-2.5 


5.2 


-107 


-3.6 


± 2.1 


-15.6 ± 2.2 


300 ± 68 


N, S 


228790063 


20:45:08.4 


-41:14:00 


-1.4 


5.0 


9 


-5.7 


± 2.2 


-6.1 ± 2.2 


170 ± 47 


S 


304920015 


21:03:42.6 


-40:23:29 


-2 


5.4 


-39 


-5.8 


± 4.7 


-10.0 ± 3.8 


152 ± 84 


S 


295010036 


21:07:15.8 


-36:06:45 


-1.2 


5.7 


-151 


2.9 


± 2.3 


-8.5 ± 2.3 


151 ± 33 


S 


304920063 


21:08:33.1 


-40:47:28 


-0.4 


6.9 


15 


4.4 


± 3.8 


-10.1 ± 3.8 


326 ± 171 


S 


295010057 


21:16:13.3 


-35:09:10 


-1.6 


5.7 


-98 


19.1 


± 2.6 


-7.6 ± 2.6 


442 ± 46 


s 


295010102 


21:25:59.7 


-36:09:49 


-1.6 


5.8 


-37 


2.1 


± 1.8 


-4.5 ± 1.8 


112 ± 47 


s 


229480002 


21:33:50.1 


-40:53:38 


-1.9 


6.2 


-56 


5.3 


± 2.5 


1.7 ± 2.5 


336 ± 81 


s 


294930009 


21:39:43.1 


-30:53:39 


-0.9 


6.4 


1 


-4.5 


± 3.1 


-6.4 ± 2.3 


168 ± 38 


s 


294930012 


21:41:20.0 


-29:53:38 


-1.6 


6.2 


-175 


9.8 


± 5.1 


-10.8 ± 5.8 


278 ± 131 


s 


295160054 


22:19:58.2 


+04:52:34 


-1.2 


7.7 


56 


-11.3 


± 2.2 


-10.1 ± 2.2 


275 ± 15 


A, T 


295160017 


22:27:58.7 


+03:59:28 


-2.2 


7.7 


-228 


-21.0 


± 5.7 


-31.0 ± 5.4 


351 ± 62 


A 


295160011 


22:28:36.2 


+06:21:09 


-2.6 


7.8 


-14 


13.6 


± 3.4 


-7.3 ± 3.4 


186 ± 29 


A, T 


303320115 


22:34:19.2 


+09:16:19 


-2 


8.4 


-73 


-8.9 


± 5.1 


11.0 ± 5.1 


463 ± 78 


A 


303320025 


22:46:15.8 


+07:51:34 


-1.5 


8.0 


-37 


0.0 


± 5.7 


-2.0 ± 5.7 


194 ± 49 


A 


303320016 


22:48:46.3 


+ 10:51:48 


-2.5 


8.7 


-180 


4.9 


± 4.8 


-25.0 ± 4.8 


386 ± 89 


A 


295130017 


23:19:16.9 


-37:40:48 


-1.6 


7.6 


94 


9.0 


± 6.3 


-9.2 ± 3.1 


134 ± 68 


S 


229410002 


23:28:24.9 


-35:56:04 


-2 


8.4 


-77 


5.7 


± 1.9 


-3.9 ± 2.0 


143 ± 55 


S 



TABLE 1. (continued) 



NAME" 


RA 
(2000.0) 


DEC 


[Fe/H] 


r b 
(kpc) 


RV C 
(km/s) 


fjL a * 

(mas/yr) 


M<5 
(mas/yr) 




V 

(km/s) 


Source d 


229410008 


23:28:40.2 


-33:55:22 


-3 


7.6 


-85 


-5.7 ± 1.9 


-11.7 ± 1 


.9 


208 ± 26 


A, S 


229410016 


23:31:27.3 


-33:57:02 





7.9 


43 


8.7 ± 3.7 


3.5 ± 3 


.2 


305 ± 70 


S 


229410022 


23:33:11.6 


-35:22:03 


-1.7 


9.0 


-25 


11.9 ± 1.8 


-2.5 ± 2 


.5 


342 ± 64 


S 


229410028 


23:35:19.6 


-36:48:25 


-2 


8.4 


-6 


6.1 ± 3.0 


-14.3 ± 3 


.0 


218 ± 89 


S 


294990022 


23:38:44.8 


-22:31:20 


-2.2 


7.8 


101 


11.2 ± 1.4 


-16.12 ± 1 


.4 


152 12 11 


A, N, S, T 


229410037 


23:38:54.2 


-35:16:38 


-2 


7.7 


-69 


2.9 ± 2.6 


-9.1 ± 2 


.9 


102 ± 35 


A, S 


294990019 


23:39:03.0 


-23:12:36 


-1.6 


9.4 


-69 


8.1 ± 4.5 


-0.4 ± 4 


.4 


254 ± 146 


S 


294990001 


23:39:19.5 


-27:24:37 


-3 


9.5 


209 


2.0 ± 2.5 


-5.4 ± 2 


.5 


238 ± 44 


S 


229660029 


23:41:15.9 


-30:24:00 


-2.5 


8.1 





17.9 ± 5.9 


-5.9 ± 4 


.3 


291 ± 111 


s 


294990035 


23:43:11.9 


-26:38:47 


-2.5 


7.9 


-203 


16.3 ± 2.6 


-19.6 ± 3 


.0 


235 ± 38 


A, S 


294990026 


23:43:31.5 


-22:57:50 


-1.6 


8.8 


-73 


10.5 ± 3.9 


-6.0 ± 3 


.7 


183 ± 101 


S 


229410054 


23:44:03.2 


-33:01:35 


-2 


9.3 


-155 


6.3 ± 4.1 


-11.1 ± 4 


.5 


259 ± 134 


S 


229660042 


23:44:29.1 


-28:19:21 


-2.1 


7.9 


-46 


4.2 ± 3.1 


-20.3 ± 2 


.5 


129 ± 44 


A, S 


228940048 


23:45:29.5 


-01:57:27 


-2 


12.1 


-135 


6 ± 5 


-4 ± 


5 


171 ± 203 


N 


294990040 


23:46:14.7 


-25:45:07 


-3 


9.3 


-66 


2.8 ± 3.5 


-8.7 ± 3 


.5 


82 ± 98 


S 


229410050 


23:46:28.1 


-34:56:01 


-1.8 


8.4 


58 


6.2 ± 5.5 


-0.6 ± 5 


.5 


207 ± 133 


s 


229660061 


23:46:48.6 


-30:00:29 


-1.8 


8.4 


-131 


-0.9 ± 7.0 


-10.5 ± 8 


.7 


181 ± 106 


s 


294990036 


23:46:51.9 


-26:48:08 


-2.7 


9.7 


48 


11.2 ± 3.3 


-9.9 ± 3 


.3 


332 ± 108 


s 


295170002 


23:47:04.3 


-16:41:01 


-1.3 


8.5 


2 


5 ± 5 


6 ± 


5 


321 ± 95 


N 


294960025 


23:47:06.2 


-30:02:49 


-3 


9.5 


8 


-1.7 ± 7.4 


-5.3 ± 8 


.0 


171 ± 264 


S 


294990042 


23:47:34.9 


-24:55:09 


-1.9 


8.5 


-29 


7.2 ± 3.6 


-4.6 ± 3 


.8 


114 ± 84 


S 


295170007 


23:47:52.7 


-14:47:59 


-2.4 


13.3 


-309 


-2 ± 5 


-7 ± 


5 


376 ± 220 


N 


294990037 


23:49:22.6 


-26:35:43 


-2.8 


8.2 


-104 


4.5 ± 3.3 


-7.5 ± 4 


.4 


104 ± 47 


S 


294990038 


23:49:40.0 


-26:31:36 


-2.4 


8.9 


-33 


13.1 ± 3.5 


-9.5 ± 3 


.4 


264 ± 95 


S 


229660059 


23:49:59.0 


-29:56:11 


-2.8 


11.5 


-126 


0.6 ± 3.2 


-5.5 ± 2 


.7 


153 ± 81 


S 


228760011 


23:51:02.7 


-34:04:00 


-1.7 


8.5 


117 


15.5 ± 4.8 


-3.5 ± 4 


.0 


343 ± 102 


S 


229660065 


23:51:16.8 


—31:30:34 


—2.9 


11.3 


238 


1.5 ± 2.9 


—4.8 ± 3 


,7 


248 ± 67 


s 


229660076 


23:51:40.0 


-29:10:16 


-0.7 


8.1 


-32 


-2.5 ± 2.3 


-7.9 ± 2 


.9 


169 ± 47 


s 


228760006 


23:52:19.2 


-36:03:54 


-2.3 


8.3 


-28 


7.4 ± 4.0 


-0.2 ± 4 


.2 


217 ± 101 


s 


228760010 


23:53:52.2 


-33:51:51 


-2.2 


9.0 


96 


6.9 ± 2.4 


-8.3 ± 2 


.5 


144 ± 54 


s 


229660069 


23:54:58.8 


-31:05:53 


-1.7 


8.3 


-114 


-8.9 ± 3.0 


-20.0 ± 1 


.9 


399 ± 54 


s 


294990064 


23:55:32.8 


-25:32:52 


-2.1 


9.7 


-48 


7.5 ± 2.2 


-4.3 ± 4 


.2 


151 ± 100 


s 


229660071 


23:55:48.3 


-30:35:43 


-3 


11.5 


35 


2.0 ± 2.8 


-6.7 ± 3 


.1 


125 ± 137 


s 


228760022 


23:58:06.4 


-33:45:11 


-2.5 


8.3 


38 


17.0 ± 2.7 


-6.7 ± 2 


.6 


278 ± 58 


s 


228760019 


23:59:20.1 


-33:17:04 


-2.1 


8.8 


2 


1.5 ± 2.9 


-12.9 ± 2 


.9 


172 ± 79 


s 



"The names of FHB stars follow Wilhelm et al. (1999). 
b heliocentric radial velocity 
c Galactocentric distance 

d A, N, S, and T denote the STARNET catalogue, the NPM catalogue, the SPM catalogue and the Tycho-2 catalogue, 
respectively. 

*the four highest- velocity FHB stars 



